Back

Proceedings of the Royal Society B: Biological Sciences

The Royal Society

Preprints posted in the last 7 days, ranked by how well they match Proceedings of the Royal Society B: Biological Sciences's content profile, based on 341 papers previously published here. The average preprint has a 0.28% match score for this journal, so anything above that is already an above-average fit.

1
Spatial variation in incidence of meningococcal meningitis: evidence from a large historical epidemic in Glasgow

Stewart, G.; Schroeder, M.; Mancy, R.; Angelopoulos, K.

2026-05-30 epidemiology 10.64898/2026.05.28.26354324 medRxiv
Top 4%
1.7%
Show abstract

Large epidemics of invasive meningococcal disease are rare in temperate regions. Here, we analyse administrative data on the largely forgotten epidemic of bacterial meningococcal meningitis that occurred in Glasgow in 1907, probably the largest on record in the UK. The epidemic, predominantly confined to the city, killed around 1,000 people, had a case fatality rate of nearly 70%, and hit infants and young children the hardest. We show the rapid rise and fall in cases and the spatial distribution of incidence and mortality rates within the city. We find that within-household overcrowding was a key driver of incidence whereas between-household geographic proximity was not. We also find that the spatial distribution of disease risk during the epidemic persisted in the post-epidemic period and during a later outbreak. The findings suggest that interventions should prioritise populations in areas that have experienced higher incidence rates to mitigate the risk of future outbreaks.

2
Fisher information matrix computation for joint longitudinal and survival models to support clinical study design and covariate effect assessment

Fayette, L.; Brendel, K.; Mentre, F.

2026-06-01 pharmacology and therapeutics 10.64898/2026.05.28.26354340 medRxiv
Top 8%
0.5%
Show abstract

Joint modelling of longitudinal data using non-linear mixed effects models and time-to-event outcomes provides a suitable framework to account for informative censoring when estimating biomarker dynamics and quantifying event risk using covariates and longitudinal trajectories. Their usefulness in clinical research depends on data collection design, particularly to precisely estimate the association (link) parameter between longitudinal and survival processes. However, optimal design strategies have so far been addressed separately for longitudinal and survival endpoints and remain unexplored for joint models. We propose two Fisher Information Matrix (FIM) computation methods for joint models, relying on Monte-Carlo integration over observations combined with either Markov Chains Monte-Carlo or Adaptive Gaussian Quadrature to integrate random effects. Their accuracy is assessed against clinical trial simulations in an oncological example based on the HORIZON III study with a tumour-growth-survival model including discrete and continuous covariates. We apply these methods to quantify the impact of follow-up duration, sampling richness, sample size, and covariate distribution on parameter uncertainty and test power. In our example, longitudinal-parameter uncertainty is barely affected by follow-up duration or sampling richness, whereas survival-parameter uncertainty decreases substantially from 1-year to 2-year follow-up. The number of subjects needed (NSN) to achieve <15\% uncertainty on the link parameter is comparable for a 2-year rich design and a 3-year sparse design. Optimal covariate distributions are stable across designs and systematically improve test power, outperforming longer and richer but non-optimised designs. These FIM-based methods accurately predict uncertainty and test powers, enabling design evaluation and NSN computation for joint-model-based clinical studies.

3
Same household, different choices: variation in health behaviors related to respiratory viruses in Illinois

Larsen, S. L.; Yang, J.; Haslett, E. M.; Anastasi, A.; Venegas, A.; Schieleit, L.; Mahmud, A.; Martinez, P. P.

2026-05-28 public and global health 10.64898/2026.05.26.26354179 medRxiv
Top 9%
0.4%
Show abstract

While SARS-CoV-2 and influenza continue to place a significant burden on population health, within-household differences in decisions towards vaccination and seeking care across these two pathogens, and across sociodemographic groups, remain largely unexplored. By conducting a household-level survey in Illinois, we found that many individuals made inconsistent decisions about vaccination: among all adults, 29% were vaccinated for only one of COVID-19 or influenza, and among those with children in the home, 39% lived with a child whose influenza or COVID-19 vaccination status differed from their own. A higher proportion of adults were vaccinated against COVID-19 compared to influenza, while the opposite was true for those younger than 18 years old. These differences hold even when accounting for disparities in coverage by age, race/ethnicity, political affiliation, and socioeconomic status. While vaccinated individuals consistently reported wanting to protect themselves or others, those who declined vaccination reported highly heterogeneous reasons ranging from resource constraints to distrust or misconceptions about vaccination. These differences are even more pronounced for COVID-19, with larger partisan gaps and higher refusal driven by safety concerns, lack of trust, or religious reasons than those who decide not to get the influenza vaccine. In contrast to vaccination, the decision to seek medical care when sick showed opposite sociodemographic trends, that are likely attributable to illness severity. Our findings highlight that closing gaps in COVID-19 and influenza vaccination coverage will require an integrative strategy that accounts for diverse motivations, fears, and barriers to access, while addressing social inequalities common to both diseases.

4
Two anti-phase spatial modes and a candidate spatial-persistence regime transition of SARS-CoV-2 in Japan: a 159-week prefecture-level sentinel surveillance study

Nakano, T.; Onozuka, D.; Ikeda, Y.; Washiyama, K.; Takashima, Y.

2026-05-26 epidemiology 10.64898/2026.05.24.26353972 medRxiv
Top 10%
0.3%
Show abstract

Background. On 8 May 2023 the Japanese Ministry of Health, Labour and Welfare reclassified COVID-19 under the Infectious Disease Control Law from a designated infectious disease (with case-by-case reporting requirements comparable to those of a Category-2 disease) to a Category-5 ("Class-5") notifiable disease, joining the same category as seasonal influenza and most other endemic respiratory infections. Under this regime, COVID-19 case counts are reported weekly from a nationwide network of sentinel medical facilities (initially approximately 5,000, reduced to approximately 3,000 following an April 2025 surveillance reform), and individual case reporting is no longer required. We aimed to characterize the spatial topology of COVID-19 epidemics under this sentinel-surveillance regime and to detect, in a data-driven manner, any structural change in epidemic dynamics over this period. Methods. We analyzed weekly per-sentinel-facility COVID-19 case counts in all 47 prefectures of Japan from 2023-W17 to 2026-W19 (159 weeks). For each week we computed the Shannon pseudo-entropy S of the prefecture-share distribution and global, local, and time-lagged Moran's I across a 92-edge contiguity-based adjacency matrix. To identify any structural change in a data-driven manner, we adopted a two-stage approach motivated by an empirical regularity established in Section 3: we first verified the wave-amplitude-invariant entropy ceiling (S_max >= 3.80 in all five pre-transition waves), then restricted change-point detection to the weeks after S(t) last attained this ceiling, applying PELT, CUSUM, and Bai-Perron sup-F within this restricted region. Seasonal structure was characterized by truncated Fourier regression with first-order autoregressive errors (Cochrane-Orcutt) over harmonic orders K = 1 to 6; between-period comparisons used moving block bootstrap as the principal inferential statistic. Results. The five epidemic waves during 2023-2025 followed a stereotyped spatial template in which S(t) traced a characteristic U-shape around each peak, with a wave-amplitude-invariant entropy ceiling reaching on average 99.4% of the theoretical maximum ln 47 (range 3.820-3.836, SD 0.006). The last week in which S(t) attained this entropy ceiling was 2025-W42. Restricting change-point detection to the 29 subsequent weeks, PELT and CUSUM localised the structural break to late 2025: PELT identified 2025-W48 (robust across penalty values >= sigma^2*ln(n) and across entropy-ceiling thresholds 3.78-3.82) and CUSUM peaked at 2025-W50 (p < 0.0001), placing the break within a two-week window centred on late November 2025. Bai-Perron sup-F peaked later at 2026-W02 (p = 0.062, with reduced power on n = 29). We adopted 2025-W48 as the principal change-point, defining 135 pre-transition weeks and 24 post-transition weeks. Two anti-phase spatial modes were identified in the pre-transition record: a summer-onset Okinawa-seeded Kyushu cascade (Mode A; annual peak epi week 26) and a winter-onset Tohoku-centred connected-cluster mode (Mode B; annual peak epi week 51), approximately 25 epi weeks out of phase. After the regime transition, this ceiling was not attained, and the spatial-persistence ratio I(tau = 8 wk)/I(0) shifted from a highly variable distribution centred near 0.27 (pre-transition, 125 weeks) to a tightly clustered distribution around 0.89 (post-transition, 24 weeks); the mean difference was 0.62 (95% bootstrap CI 0.32 to 0.90; moving block bootstrap p < 0.0001 across block lengths 1-12). The principal finding remained significant under autoregressive-augmented null models and was robust to adjacency-matrix choice, the April 2025 surveillance reform, harmonic order K = 1 to 6, and Okinawa exclusion. Conclusions. Data-driven analysis of 159 weeks of Japanese sentinel surveillance identifies a candidate spatial-persistence regime transition emerging in late November 2025, in which the spatial structure of weekly case shares persists for at least 8 weeks rather than dissipating as in pre-transition. The transition coincides with loss of the wave-amplitude-invariant entropy ceiling and with absence of the Mode A signature through the observed post-transition period. The recent uptick in Okinawa case shares (continuing through 2026-W19) leaves open whether the Mode A signature is structurally suppressed or merely deferred; observation through summer 2026 is required to distinguish a sustained shift from a transient anomaly.

5
Keeping human in the loop: A three-phase generative AI workflow for research integrity in data-intensive science.A methodological case study using elite Ethiopian distance-running data

Galko, P.; Yisamaw, A.; Haugen, T.; Seiler, S.

2026-05-29 sports medicine 10.64898/2026.05.29.26354013 medRxiv
Top 10%
0.3%
Show abstract

Background: Generative AI tools can support data-intensive research by writing code, drafting prose, searching analytical possibilities, and stress-testing claims. They can also produce false citations, drift between statistical specifications, and lose continuity across long investigations. This paper describes a practical workflow for using AI systems in empirical research while keeping discovery, verification, and accountability inspectable. Methods: We developed and applied a three-phase human-AI workflow to a case study of 14 elite Ethiopian distance runners. The dataset contained 22,605 GPS-segments collected across 97 consecutive days in late 2025, supplemented by venue and athlete metadata collected in the field. Phase 1 used an autonomous data-exploration tool to pre-filter the hypothesis space across five seeded research questions. Phase 2 used an AI system under direct human guidance to construct candidate findings into numerical claims, verification scripts, and draft text. Phase 3 used an independent AI system in an adversarial role to stress-test methods, statistics, prose, figures, and citations. The workflow was informed by Pearl's distinction between association, intervention, and counterfactual reasoning, with human judgement retained for research direction, interpretation, and final claims. Results: The workflow produced three empirical analyses and a documented correction process. The analyses estimated an altitude-to-sea-level pace correction of +0.10 min/km per 1,000 m at matched heart rate, showed why pooled altitude-surface regression was not identifiable within this venue system, documented method-dependence in heart-rate-based intensity classification, characterised within-venue route variation as a 64/36 path-fixed-to-trail-variable split with the Sululta label resolving into two functionally distinct sub-venues, and reframed the cohort's training through a 3x3x3 prescription lattice grounded in Ethiopian coaching practice. The adversarial phase identified several hallucinated citations, a terminology error between HC1 and cluster-robust standard errors, and several inconsistencies between prose, figures, and computed results. Verification scripts re-derived nearly all numerical claims from the cleaned lap-level data. Conclusions: The case study shows how researchers can organise AI-assisted empirical work so that candidate discovery, claim construction, independent stress-testing, and final accountability remain separated. The workflow did not remove the need for domain expertise or human judgement. Its value was in making the route from candidate finding to manuscript claim explicit, reproducible, and open to challenge. Trial registration: Not applicable.

6
Intention of UK residents to wear facemasks and practise social distancing during the next respiratory virus pandemic

Smith, D. R.; Buckell, J.; Hancock, T. O.; Morrell, L.; Pouwels, K.

2026-05-30 public and global health 10.64898/2026.05.21.26353824 medRxiv
Top 12%
0.3%
Show abstract

Background: Wearing facemasks and practising social distancing slow the spread of respiratory pathogens. However, in the event of a new pandemic emerging, the willingness of populations to voluntarily adopt these behaviours is unclear. Methods: A discrete choice experiment was conducted among 2,006 UK-based adults. Participants were presented with hypothetical scenarios describing the emergence of a respiratory virus pandemic and were asked to choose when they would wear facemasks and practise social distancing. A mixed multinomial logit model was used to jointly estimate how disease severity and prevalence, uncertainty in these quantities, and individual-level characteristics influence behavioural choices. Findings: Participants were averse to facemasks and social distancing in the absence of pandemic risk. For each ten-unit increase in severity (10 additional hospitalisations/1,000 infections), the odds of always wearing a facemask outside the home increased by 15.9% (95%CI: 14.3%, 17.5%), relative to rarely/never, and the odds of avoiding all people as much as possible increased by 16.4% (14.6%, 18.2%), relative to not avoiding anyone. Greater disease prevalence, uncertainty in disease severity or disease prevalence, a university education, prior COVID-19 vaccination and non-white ethnicity were also associated with choosing to always wear facemasks and avoid all people as much as possible. The probability of participants choosing to rarely/never wear facemasks varied from 13.4% (11.9%, 14.9%) in the lowest-risk scenario to 1.4% (1.2%, 1.7%) in the highest-risk scenario. Interpretation: Perceived risks of disease and associated uncertainty drive intention of UK adults to adapt their behaviour in a future pandemic.

7
Why epidemic risk at the 2026 World Cup may not be what you think

Lessler, J.; Smith, C. P.; Das, P.; Sykes, A. L.; Urbinati, A.; Geith, K.; Powers, K. A.; Davis, J. T.; Kern-Allely, S. C.; Vega Yon, G. G.; Lofgren, E. T.; Pearson, C. A. B.; Vespignani, A.

2026-06-01 epidemiology 10.64898/2026.05.28.26354384 medRxiv
Top 13%
0.2%
Show abstract

Background: The 2026 FIFA World Cup may bring over one million visitors to North America from around the globe to participate in mass gathering events. The nature of the event and recent news have raised concerns for some that the tournament could lead to infectious disease outbreaks or fuel existing epidemics. Objective: To systematically assess the infectious disease threat posed to the United States by the tournament. Design: A multi-institutional team evaluated pathogen-specific risk across three dimensions: importation, outbreak potential, and impact to identify a priority pathogen list. A systematic screening protocol ensured common criteria and that pathogen information was collected when necessary to inform inclusion. Results: Increased risk from the World Cup is near zero for 63 of 77 evaluated pathogens. Pathogens were predominantly excluded as threats due to low excess importation risk and low outbreak potential if introduced. The remaining priority pathogens fall into five categories: (a) mosquito borne pathogens with the potential for sustained transmission in some host cities, (b) seasonal respiratory viruses, (c) chronic infections with high prevalence outside the United States, (d) pathogens present in the United States with likely increased transmission at World Cup activities, and (e) high-consequence infectious threats. Limitations: Data availability is variable across diseases. Impact calculations may not reflect actual costs to host cities. Disease incidence in World Cup travelers may differ from national incidence rates. Conclusion: While infectious disease outbreaks at the 2026 FIFA World Cup are possible, in an already highly connected world where large gatherings are frequent, the elevated risk from the tournament is not as extreme as it first may seem.

8
Inferring Sexual Network Bridging Using Genomics: A Simulation Study

Kline, M. C.; Helekal, D.; Oliveira Roster, K. I.; Grad, Y.

2026-05-26 infectious diseases 10.64898/2026.05.24.26353967 medRxiv
Top 13%
0.2%
Show abstract

The dynamics of sexually transmitted infections involve interconnected transmission networks, including men who have sex with men and heterosexual populations. Understanding the extent of bridging between these networks can inform surveillance, guide interventions, and aid in the interpretation of their impact, but methods for quantifying bridging have been lacking. Here, we addressed whether pathogen genomics tools, successfully used to reconstruct transmission in other contexts, could accurately infer sexual network bridging. Based on simulations of gonorrhea spread, we evaluated phylodynamic bridging metrics inferred by ancestral state reconstruction under a range of sampling schemes, from comprehensive to sparse. These metrics differentiated sexual network structures even with biased sampling schemes, but accuracy depended on the sampling scheme and density: phylodynamic bridging estimates using sequences from all detected infections for one network configuration were on average 6.9% above the true value, whereas estimates from 5% of infections in symptomatic men with many partners were on average >1000% above the true value. These results suggest routine overestimation of bridging from unadjusted inferences from genomics data and provide context for interpreting existing genomic surveillance data and targeted studies.

9
Data Assimilation Substitutes for Biological Complexity in Hybrid Influenza Forecasting Models

Alleman, T. W.; Van Wesemael, T.; Shanker, N.; Mietchen, M. S.; Loo, S.; Ajagbe, S. O.; Baetens, J. M.; Lemaitre, J.; Hill, A. L.; Truelove, S. A.; Bento, A. I.

2026-05-27 public and global health 10.64898/2026.05.19.26353597 medRxiv
Top 14%
0.2%
Show abstract

Hybrid mechanistic-statistical models offer interpretability and adaptability for short-term seasonal epidemic forecasting, but it remains unclear whether their accuracy depends more on increased biological complexity or on the assimilation of richer data. Using eight retrospective influenza seasons in North Carolina, we evaluate whether training on historical data and assimilating auxiliary emergency department (ED) visit data improves four-week-ahead hospital admission forecasts more than adding biological complexity (multi-subtype structure and cross-season immunity). Hierarchical Bayesian training on historical data improves accuracy by 22.4 % (95 % CI: 16.4-28.1 %), and inclusion of ED visit data yields a further 5.3 % (95 % CI: 3.0-7.6 %) improvement, whereas added biological complexity produces diminishing or null gains. We further observe a substitution effect in which ED visit data partially compensates for omitted biological structure. We deployed a simplified model variant in the 2025-2026 CDC FluSight Challenge and ranked among the top ensemble performers, supporting the robustness of Bayesian hierarchical training in real time. Together, these findings indicate that short-term forecast accuracy is driven more by historical learning and assimilating auxiliary signals than by biological fidelity, with implications for how forecasting systems should balance mechanistic complexity.

10
HIV-HPV Syndemic and Anal Precancerous Lesions Among MSM and Transgender Women in Pakistan: A Biological Continuum in High-Risk Sexual Networks

Ejaz, M.; Ahmed, A.; Rizvi, S. H.; Rizvi, A. A.; Ali, F.; Haroon, A.

2026-06-01 public and global health 10.64898/2026.05.28.26354356 medRxiv
Top 15%
0.2%
Show abstract

Background: Sexual and gender minorities (SGM), including men who have sex with men (MSM) and transgender women, often face stigma, legal constraints, and limited access to sexual and reproductive health services. These conditions restrict prevention and care, increasing vulnerability to HIV and human papillomavirus (HPV) infections. While strong HIV-HPV interaction is documented in high-income settings, evidence from low- and middle-income countries remains limited. This study examines the burden, co-infection dynamics, and progression of HPV infection and anal dysplasia among MSM and transgender women in Pakistan. Methods: A cross-sectional study was conducted between September 2015 and October 2016 among men who have sex with men (MSM) and transgender women recruited from sexual health and antiretroviral therapy centers in Karachi. Eligible participants were aged [&ge;]18 years and self-reported anal sex within the past 6 months (N=298). Two anal specimens were collected for HPV DNA detection and genotyping using PCR, and anal squamous intraepithelial lesions (ASIL) were assessed cytologically using the Bethesda classification. Associations were estimated using Cox proportional hazards regression algorithms to derive prevalence ratios (PRs). Results: Among participants, 44% (n=133) were living with HIV. Overall HPV prevalence was 65.1%, rising to 87% among HIV-positive individuals compared to 48% among those without HIV ({chi}{superscript 2}p[&le;]0.001). Likewise 28.9% of participants living with HIV were infected with two or more than two types of HPV as compared with 18.8% participants without HIV ({chi}{superscript 2}p[&le;]0.001). HIV infection was strongly associated with HPV acquisition (adjusted PR 2.81, 95% CI 2.16-3.82). Among HPV-positive participants (n=194), 58.8% were co-infected with HIV. High-risk HPV was highly prevalent among those living with HIV (83.2% vs. 35.3% ({chi}{superscript 2}p[&le;]0.001)), with HPV16 as the dominant oncogenic type. Multiple HPV infections were more common among HIV-positive individuals ({chi}{superscript 2}p[&le;]0.001), and HIV seropositivity was 3.43 (95% CI: 2.55-3.51) times higher among those with high-risk HPV. Co-infected participants demonstrated prolonged smoking, longer duration of sex work, high-intensity sex work with limited condom negotiation, and higher prevalence of anal warts (all p<0.05). Anal dysplasia (ASIL) was present in 35% of participants and was higher among HIV-positive individuals (42.4% vs. 28.1%, p<0.001). HIV-HPV co-infection was independently associated with ASIL (adjusted PR 1.75, 95% CI 1.07-2.88), while high-risk HPV further amplified this risk (PR 3.04, 95% CI 1.75-5.26). Conclusion: These findings demonstrate a biological continuum in HIV-positive MSM and transgender women, where HIV increases HPV acquisition, persistence, and multiplicity, accelerating progression to anal dysplasia. This reflects a syndemic shaped by biological interaction and structural vulnerability. Integrating HPV screening and vaccination within HIV services is essential to interrupt progression to cancer in this high-risk population.

11
High-dimensional Characterization of Genome-Environment Fitness Landscapes in Klebsiella pneumoniae

Zhou, G.; Williams, G.; Millner, M. T.; AlHirayban, R.; Alosaimi, W.; Fallatah, O.; Hart, A. J.; Malaikah, M.; Iftikhar, S.; Ahmad, H.; Roghanian, M.; Mustonen, V.; AlYami, R.; Banzhaf, M.; Moradigaravand, D.

2026-05-30 genetic and genomic medicine 10.64898/2026.05.28.26354339 medRxiv
Top 15%
0.2%
Show abstract

Background Bacterial fitness is shaped by interactions between genome variation and environmental context, yet how these interactions determine its predictability and heritability remains unclear. In the clinically important pathogens of Klebsiella pneumoniae, a leading cause of hospital-acquired infections, this question is particularly pressing. Despite extensive genomic characterization, we still lack a systematic understanding of how genome-wide variation translates into fitness across diverse environments in K. pneumoniae. Methods We filled this gap by profiling a systematic collection of 1,462 clinical K. pneumoniae isolates across 214 diverse environmental and pharmacological stress conditions using high-throughput chemical genomics. Fitness was quantified from colony growth and integrated with whole-genome sequencing data. Genome-wide association analyses identified genetic determinants of fitness, and machine learning models incorporating genomic features were used to predict fitness.Results Fitness exhibited a strongly environment-dependent genetic architecture, with modest but significant concordance between genetic background and phenotypic variation. Under antibiotic and stress-combination conditions, fitness was driven by discrete, high-effect determinants, including known resistance genes, resulting in stronger signals and improved predictability. In contrast, non-antibiotic environments showed more polygenic and distributed architectures with weaker associations. Genome-wide analyses identified both established and previously uncharacterized genes linked with fitness across conditions. Resistance and virulence determinants exhibited clear context-dependent trade-offs, conferring fitness advantages under selection but imposing costs in non-selective environments. Consistent with this, plasmid carriage showed environment- and genotype-dependent fitness effects, with benefits under antibiotic pressure and measurable costs otherwise. Genomic variant-based models for fitness prediction achieved moderate performance (Mean Spearman correlation ({rho}) = 0.36 (95% CI: 0.18-0.67) for predicted versus observed values in unseen data) across conditions, with improved accuracy under strong antibiotic selective pressures, and produced well-calibrated prediction intervals with high coverage. Despite strong population structure effect on predictions, models captured predictive gene and SNP biomarkers for fitness. Conclusion These findings highlight that bacterial fitness is an emergent property of genome-environment interactions rather than a fixed attribute of genotype. This work establishes a unified high-dimensional genotype-phenotype framework linking genomic variation to fitness across diverse conditions in a major pathogen, with broader implications for other pathogenic bacterial species.

12
Consumer Opinions, Lot-to-Lot Variability, and Pharmacokinetics of Transdermal Melatonin Products: A Randomized, Crossover Clinical Trial

Bonilla, K.; Sherman, V. M.; Arbaiza, A. S.; Dougherty, M.; Olson, L. E.

2026-05-29 pharmacology and therapeutics 10.64898/2026.05.27.26354234 medRxiv
Top 16%
0.1%
Show abstract

In some countries, melatonin is sold without a physician prescription and dosage is unregulated. Transdermal products have become popular including those marketed for children. We measured consumer assumptions about these products among adult residents of the United States, analyzed lot-to-lot variability, and compared the pharmacokinetics of melatonin administered in oral, lotion, and bath product forms. Survey respondents (n=199) believed oral melatonin was more effective than transdermal products and that all melatonin products were relatively safe. Melatonin lotion products analyzed by HPLC displayed lot-to-lot variability as well as changes in formulation and product claims. To determine pharmacokinetics, three different treatments (oral tablets, lotion, and bath immersion) were administered to twelve undergraduate participants in a randomized, crossover design. Five additional participants completed bath product treatment only. Participants collected saliva samples up to 48 hours after administration, which were analyzed for melatonin by enzyme-linked immunosorbent assay. Oral (n=11) and lotion formulations (n=12) caused maximum salivary melatonin levels within 30 minutes after administration, but bath immersion did not cause increases in saliva melatonin (n=17). The half-life of oral melatonin was 1.17 [0.69 -- 1.65] hours versus 5.72 [3.75 -- 7.68] hours for lotion treatment (p = 0.011, effect size r = 0.770). Melatonin lotion may pose a risk to consumers who assume it is safe and less effective than oral tablets, when in fact it may be very potent and remain at high physiological levels into the following day. This study is registered on clinicaltrials.gov (NCT06382610) and was funded by the Sleep Research Society.

13
Cross-Sectional Measures of Periodontal Severity: Distortion from Severity-Dependent Tooth Loss

McCormick, K. M.; Amarasena, N.; Guzzo, G.; Nath, S.; Jamieson, L.

2026-05-30 dentistry and oral medicine 10.64898/2026.05.27.26354277 medRxiv
Top 16%
0.1%
Show abstract

Aim: Cross-sectional summaries of periodontitis based on clinical attachment loss (CAL) are, by definition, conditioned on surviving teeth. Because the most severely affected teeth are more likely to have been lost, these measures may underestimate cumulative disease burden and show an artificial flattening (attenuation) of severity with age. We hypothesised that measures more sensitive to severe attachment loss would show greater attenuation at older ages than measures defined across a broader range of sites. Materials and Methods: Using nationally representative data from adults aged 30+ years in NHANES 2009-2014, we examined age-specific trajectories across multiple continuous measures of periodontal severity and assessed whether divergence between measures followed the pattern predicted under severity-dependent tooth loss. Results: The proportion of observable sites declined from 93% at ages 30-34 to 68% at 80+ years, establishing the structural basis for the divergence observed across severity measures. All severity measures showed nonlinear attenuation with age, with distortion increasing with severity threshold. Higher-threshold measures exhibited the greatest attenuation, while lower-threshold measures showed more stable trajectories. Conclusions: Cross-sectional summaries of periodontitis reflect disease among surviving teeth rather than cumulative damage across teeth originally at risk. Attenuation at older ages is consistent with depletion of the most severely affected teeth rather than biological slowing. Distortion varies by measure, with higher-threshold and mean-based indices most affected, whereas the CAL 3+ mm threshold provides a more stable basis for age comparisons.

14
Widespread Self-Medication and Unsafe Access to Analgesics and NSAIDs in Urban Conakry, Guinea: Prevalence, Associated Factors, Risk Exposure Profiles, and Health-System Implications. A Cross-Sectional Study of 1,032 Participants.

LAWA GARANDJI, D.; BALDE, A. O.

2026-05-30 pharmacology and therapeutics 10.64898/2026.05.21.26353180 medRxiv
Top 16%
0.1%
Show abstract

ABSTRACT Background: Self medication with analgesics and non steroidal anti inflammatory drugs (NSAIDs) is common in low- and middle income countries and may expose users to preventable adverse outcomes. Evidence from Guinea remains scarce. This study aimed to estimate the prevalence of self medication with analgesics and NSAIDs among pharmacy clients in urban Conakry, identify associated factors, and describe clinical risk situations. Methods: We conducted a pharmacy based analytical cross sectional study in 30 private pharmacies across Conakry, Guinea. A total of 1,032 participants seeking analgesics or NSAIDs were enrolled between November 3, 2012, and April 5, 2013. Self-medication was defined as acquisition or use without a valid medical prescription. Factors associated with self-medication were analysed using multivariable logistic regression. Results: Among 1,032 participants, 603 reported self medication (prevalence 58.4%). Previous unsupervised use was reported by 78.7%. The most frequently used medicines were paracetamol (56.9%, n=587), diclofenac (21.3%, n=220), ibuprofen (17.9%, n=185), and aspirin (3.9%, n=40). Overall, 68.0% (n=702) reported no knowledge of potential adverse effects. Clinical risk situations were frequent: gastrointestinal disorders (41.3%, n=426), hypertension (9.2%, n=95), and pregnancy exposure among reproductive age women (26.0%). In multivariable analysis, self medication was independently associated with previous analgesic/NSAID use (aOR = 2.8, 95% CI: 2.1 to 3.6), lack of knowledge of adverse effects (aOR = 1.9, 95% CI: 1.4 to 2.5), informal occupation (aOR = 1.6, 95% CI: 1.2 to 2.2), and age 18 to 59 years (aOR = 1.5, 95% CI: 1.1 to 2.1). Conclusions: In this pharmacy based study conducted in urban Conakry, self medication with analgesics and NSAIDs was common and frequently associated with limited awareness of potential adverse effects. These findings support the need for strengthened pharmaceutical regulation, pharmacist-led counselling, health literacy interventions, and improved access to primary care. Keywords: self medication; analgesics; NSAIDs; paracetamol; diclofenac; ibuprofen; pharmacy; Guinea; Conakry; drug safety; public health.

15
Dentine markers of pre/early postnatal lead exposure links with brain, cognitive, and behavioral outcomes in adolescents

Marshall, A. T.; Kan, E.; Adise, S.; König, M.; McConnell, R.; Martinez, M.; Midya, V.; Arora, M.; Sowell, E. R.

2026-05-27 pediatrics 10.64898/2026.05.26.26354134 medRxiv
Top 17%
0.1%
Show abstract

Lead is a toxic metal ubiquitous in our environment. While dramatic reductions in lead sources have paralleled equivalent decreases in lead-poisoning rates, chronic lead exposure remains a critical public health concern. Childhood lead exposure (at its lowest levels) is liked to changes in cognitive development but less is known about lead's effects on children's brain structure, especially as a result of in utero exposure. We measured prenatal and early-postnatal lead exposure in shed deciduous teeth of 448 9- and 10-year-old children (from 20 United States cities) and linked those lead levels to childhood brain structure, cognition/behavior, and neighborhood- and family-level socioeconomic characteristics. Here we show negative associations between tooth-lead levels and the thickness of the brain's cortex, particularly in regions linked to language processing. With increasing tooth-lead levels, children of lower-income (versus higher-income) families showed steeper declines in receptive vocabulary. Caregiver-reported behavioral problems exhibited similar associations. With in utero exposure linked to adverse neurodevelopmental outcomes (well before lead exposure and its risks are evaluated by healthcare professionals), prenatal screening of maternal lead levels/exposure, coupled with recommended strategies to reduce its placental transmission, may help reduce lead's effects on future generations.

16
Estimating Lifetime Periodontal Burden Under Informative Tooth Loss

McCormick, K. M.; Amarasena, N.; Guzzo, G.

2026-05-30 dentistry and oral medicine 10.64898/2026.05.27.26354300 medRxiv
Top 17%
0.1%
Show abstract

Background: Periodontitis is defined by cumulative, irreversible tissue destruction, yet population-based measurement typically relies on cross-sectional indicators derived from retained teeth. Destruction that occurred earlier in life, particularly disease severe enough to result in tooth loss, is structurally excluded from these measures, potentially leading to systematic underestimation of lifetime periodontal burden. Objective: To develop and evaluate a measurement framework that estimates lifetime periodontal burden from cross-sectional data by explicitly incorporating informative tooth loss under etiological uncertainty. Methods: Data were drawn from 10,324 adults aged [&ge;]30 years participating in the 20090-2016 National Health and Nutrition Examination Survey (NHANES) who completed full-mouth periodontal examination and glycated hemoglobin (HbA1c) testing. Lifetime periodontal burden was estimated by combining observed clinical attachment loss in retained teeth with probabilistic contributions from missing teeth, using three alternative age-stratified attribution schedules derived from epidemiological studies of periodontal extraction. Performance was compared with conventional measures of periodontal severity and extent using distributional analyses, correlations with HbA1c, discrimination of diabetes status, and relative importance analysis. Age-adjusted models were treated as sensitivity analyses. Results: Estimated lifetime periodontal burden exhibited strong, monotonic age gradients across glycemic categories, in contrast to more attenuated patterns observed for severity and extent. Across attribution schedules, lifetime burden showed stronger correlations with HbA1c ({rho} = 0.30-0.32) than conventional measures. In multivariable models including all indices, lifetime burden retained an independent association with HbA1c, whereas severity and extent contributed little unique information. Discriminative performance for diabetes status was consistently higher for lifetime burden than for conventional measures and remained stable across attribution schedules. Conclusions: Lifetime periodontal burden can be estimated from cross-sectional data by explicitly modelling informative tooth loss rather than restricting measurement to retained teeth. Incorporating historical tissue loss under uncertainty yields a more coherent representation of cumulative periodontal destruction than snapshot-based measures and provides a methodological basis for life-course-oriented periodontal epidemiology.

17
Using Bayesian Evidence Synthesis to estimate the number of sex workers in the United Kingdom

Long, H.; Gada, L.; Murray, L.; Laurence, T.; Hayward, A.; Finnie, T.

2026-05-26 public and global health 10.64898/2026.05.21.26353767 medRxiv
Top 17%
0.1%
Show abstract

Sex work is diverse and includes a broad range of people and settings. Over the last thirty years, a large proportion of public health emergencies of international concern (PHEIC) have involved infections transmitted through sexual or close contact and in sexual networks (WHO 2024). Sex workers can face increased disadvantage in relation to these public health emergencies. Given the significant health inequalities sex workers can face, they should be eligible to receive targeted and tailored health support to reduce health protection risks (Hester 2019; Jeal and Salisbury 2004a). However, they are often not explicitly eligible for targeted and tailored support due to a lack of information on incidence, prevalence of disease, and even more basic data such as reliable estimates of the number of sex workers in the UK. Accordingly, the aim of this paper is to determine a population size estimate, with uncertainty, that is more robust than those currently available. In this study, we apply Bayesian Evidence Synthesis to bring together historic estimation efforts with recent ONS National Population Estimates and Genito-Urinary Medicine Clinics Attendance Data (GUMCAD) from the UK Health Security Agency (UKHSA). A key feature of our model is the embedding of uncertainty from each input study in model priors, hence propagating it through to our final estimate. The Bayesian evidence synthesis model estimated a total of 84,000 sex workers in the United Kingdom (95% credible interval: 49,000-130,000), representing 0.121% of the current UK population.

18
Coaching for quality improvement under performance-based contracting: a theory-of-change evaluation in Honduras

Munar, W. J.; Aranda, L. E.; Lauria, M. E.; Bernal Lara, P.; Innocenti, C.; Rodriguez, M.

2026-05-30 health systems and quality improvement 10.64898/2026.05.21.26353487 medRxiv
Top 17%
0.1%
Show abstract

Introduction. Practice coaching is increasingly used to strengthen quality improvement (QI) capacity in primary healthcare (PHC) systems in low and middle income countries (LMICs), yet the causal pathways through which it shifts provider behaviour, and the systemic conditions that enable or constrain those pathways, remain under theorised. Using a theory based qualitative evaluation, we examined how and why a practice coaching intervention influenced QI in cervical cancer screening (CCS) and antenatal care (ANC) within Honduras decentralised PHC system during the third phase of the Salud Mesoamerica Initiative (SMI). Methods. We conducted a within case explanatory case study. A programme theory was reconstructed before data collection and iteratively refined against evidence. Data comprised semi structured interviews with 11 midlevel managers, 6 PHC team medical leads, and 2 regional managers, complemented by direct observation and document review. We applied combined deductive and inductive coding, thematic analysis, and pattern matching, and reporting per COREQ. Results. We identified four causal patterns that refined the initial programme theory. Three were activated pathways: (1) novel professional identity among participating managers; (2) collective efficacy and data driven learning, sustained through verifiable progress on observable indicators, strong for CCS but null for ANC, where outcomes were less attributable to teams actions; and (3) relational coordination, psychological safety, and trust, which provided the interpersonal basis for the first two. A fourth, unanticipated pattern showed structural misalignment between coaching enabling, learning based logic and the directive, punitive logic of Honduras performance based contracting environment, confining gains to localised enabling bubbles. Conclusion. Coaching can activate meaningful QI pathways in LMIC primary care, but sustained, equitable impact requires deliberate alignment between coaching learning oriented principles and the institutional performance management architecture, and matching of coaching investment to clinical processes with observable, attributable outcomes.

19
Sleep Disorders Modify the Age-Related Trajectory of Circadian Rest-Activity Rhythms: Evidence from NHANES 2011--2012 Wrist Actigraphy

Yin, L.; Lee, C. W.; Wong, A.

2026-06-01 epidemiology 10.64898/2026.05.28.26354369 medRxiv
Top 17%
0.1%
Show abstract

Background: Circadian rest-activity rhythms weaken with age, but whether sleep disorders modify this trajectory is unknown. Methods: We analyzed wrist accelerometry data from 4,386 participants aged 6-80 years in the 2011-2012 National Health and Nutrition Examination Survey (NHANES). Circadian features were extracted using cosinor analysis and nonparametric methods; a Circadian Disruption Index (CDI) was constructed from five standardized components. Survey-weighted regression with natural cubic splines and Wald F-tests tested age-by-sleep-disorder interactions using Taylor series linearization for variance estimation. Results: Doctor-diagnosed sleep disorder (N = 360, 8.2%) was associated with significantly different age-related trajectories of amplitude (F(2,17) = 11.24, p = 0.0008) and MESOR (F(2,17) = 8.22, p = 0.0032), both surviving Bonferroni correction (p < 0.006). CDI was higher in those with a sleep disorder (0.290 vs. 0.131, p < 0.001) and was independently associated with higher BMI (beta = 1.33 kg/m2, p < 0.001), higher HbA1c (beta = 0.089%, p = 0.004), greater diabetes prevalence (beta = 3.8 percentage points, p < 0.001), and worse depressive symptoms (beta = 0.43 PHQ-9 points, p = 0.020). Sensitivity analyses using a broader sleep problem exposure did not replicate these interactions. Conclusions: Doctor-diagnosed sleep disorders are associated with an altered age-related decline in circadian amplitude and mean activity level. CDI was independently linked to cardiometabolic and depressive outcomes, supporting a mechanistic connection between clinically significant sleep pathology and circadian disruption across the lifespan.

20
Impact of the Management Development Programme (MDP) on primary health care manager competencies and organisational Performance

Sineke, T.; Shumba, K.; Moolla, A.; Mongwenyana-Makhutle, C.; Hongoro, D.; Miot, J.; Kruger, P.; Graven, J.; Onoya, D.

2026-06-01 health systems and quality improvement 10.64898/2026.05.28.26354357 medRxiv
Top 20%
0.1%
Show abstract

Primary healthcare (PHC) managers are central to the functioning of South Africas healthcare system, yet many assume leadership roles without formal management training. To address this gap, the Aurum Institute developed the Management Development Programme (MDP), a structured leadership and management training intervention aimed at strengthening PHC management competencies. This study evaluated the impact of the MDP on leadership practices, organisational readiness for change, and workplace stress among PHC managers in the Western Cape Province. A non-randomised matched cluster trial was conducted across 20 PHC facilities. Intervention facilities were purposively selected based on participation in the MDP, while matched control facilities were randomly selected. Data were collected using structured and semi-structured surveys administered to facility managers and clinic staff. Leadership competency was assessed using the Leadership Practices Inventory (LPI), which measures five dimensions of exemplary leadership: Model the Way, Inspire a Shared Vision, Challenge the Process, Enable Others to Act, and Encourage the Heart. Organisational readiness for change was measured using Kotters 8-Step Framework, while workplace stress was assessed using a 13-item version of the Brief Job Stress Questionnaire focusing on Job Meaning, Environmental Quality, Autonomy, and Control. Intervention effects were estimated using generalised linear models adjusted for manager age, years in role, matched-pair fixed effects, and cluster-robust standard errors. Outcomes were reported as adjusted risk differences with 95% confidence intervals and two-sided p-values. A total of 20 facility managers (median age 51 years; IQR 42-55; 90% female) and 105 clinic staff members (median age 42 years; IQR 35-50) participated in the study. Managers in both intervention and control facilities reported consistently high self-rated leadership competency scores across all LPI domains, with no statistically significant differences between groups. Similarly, clinic staff rated managers highly across the standard LPI domains, and no significant differences were observed between intervention and control facilities. Despite the absence of significant differences in overall leadership competency scores, staff in intervention facilities reported significantly stronger relational and communication practices among managers compared with staff in control facilities (72.7% vs. 64.0%; adjusted risk difference 22.0%, 95% CI 6.1-37.8; p=.007). After adjustment for age and tenure imbalances, intervention facilities also demonstrated significantly higher scores for institutionalised capability and learning culture (adjusted risk difference 21.3%, 95% CI 0.6-42.0; p=.043). Managers who participated in the MDP further reported stronger perceptions of district support, including improved internal leadership and cultural readiness (adjusted risk difference 22.1%, 95% CI 14.0-30.3; p<.001) and greater district leadership and resource availability (adjusted risk difference 28.1%, 95% CI 15.6-40.6; p<.001). No statistically significant differences were observed in workplace stress across any domain. Although the MDP did not produce measurable short-term improvements in managers self-rated leadership competencies or standard LPI domains as assessed by staff, it was associated with important gains in relational leadership practices, organisational readiness for change, and perceived district support. These findings suggest that structured management training programmes may strengthen critical organisational and interpersonal foundations necessary for sustained performance improvement within PHC settings.